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agreement  with  the  asymptotic  results.  „ 
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BREAKING  SOLUTIONS  OF  THE  VARIABLE  COEFFICIENT 
KORTEWEG-deVRIES  EQUATION 

1 .  Introduction 

One  of  the  most  important  problems  in  physical  oceanography  is  the  description 
of  surface  gravity  waves  as  they  propagate  in  water  of  variable  depth.  An 
extremely  useful  tool  in  the  study  of  this  phenomenon  is  the  variable  coefficient 
Korteweg-deVries  equation.  This  equation  contains  the  most  important  effects 
of  modulation  due  to  the  depth  variation,  nonlinearity  and  dispersion.  It  has 
been  shown  that  this  equation  can  describe  the  gradual  modulation  of  solitary 
waves  and  cnoidal  wavetrains,  see  Johnson  (1973b),  Ostrovskiy  and  Pelinovskiy 
(1970)  and  Svendsen  and  Hansen  (1978),  as  well  as  the  shoal-induced  fissioning 
of  solitary  waves,  see,  e.g.,  Johnson  (1972).  In  the  appropriate  long 
wave  limit,  it  can  also  describe  waves  which  steepen  and  break.  The  main 
objective  of  this  report  is  to  describe  another  class  of  solutions  to  this 
equation.  Over  most  of  their  propagation,  these  solutions  may  be  described 
by  the  linearized  version  of  this  equation.  Near  the  shore,  the  linear  theory 
breaks  down  and  in  this  region  nonlinear  effects  predominate  and  the  wave  is 
observed  to  break  by  plunging.  The  theory  presented  here  gives  the  breaking 
point  and  the  wave  amplitude  at  breaking  explicitly  in  terms  of  the  initial, 
i.e.,  off-shore,  wave  shape  and  the  bottom  topography.  Once  the  wave  breaks  and 
a  bore  forms,  the  resultant  interaction  of  the  bore  with  the  flow  behind  it 
will  tend  to  decrease  the  bore  amplitude;  the  results  obtained  here  are  also 
expected  to  give  reasonable  order  of  magnitude  estimates  for  the  maximum 
observed  waveheight. 

Although  the  starting  point  for  our  analysis  is  the  Korteweg-deVries  equation, 
the  results  obtained  will  also  be  applicable  to  actual  flow  problems.  The 
conditions  under  which  this  is  true,  i.e.,  under  which  the  Korteweg-deVries 

equation  is  valid,  are  discussed  in  Section  2. 
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An  interesting  result  which  comes  out  of  our  analysis  is  the  dependence 

of  the  wave  shape  and  break  point  on  the  initial  parameters  of  the  problem. 

For  the  cases  considered  here,  these  depend  on  the  wave  shape  and  bottom 

€  _  u2 

variation  through  €,  u  and  5  and,  in  particular,  the  ratios  e  =  and  v  =  -r—  ; 
see  Section  3  for  this  notation.  The  results  of  numerical  or  experimental  studies 
can  therefore  be  presented  in  terms  of  these  natural  parameters. 

It  is  important  to  note  how  the  present  investigation  differs  from 
previous  studies  of  shoaling  water  waves.  While  the  authors  mentioned  above 
calculate  the  behavior  of  a  single  type  of  wave  as  it  approaches  a  beach,  we 
have  shown  how  a  wave  transforms  from  a  linear  dispersive  wave  into  one  which 
is  basically  nonlinear  and  nondispersive.  Furthermore,  we  do  not  claim  to 
solve  the  very  difficult  problems  of  run-up  and  backwash  discussed  by  many 
authors,  see,  e.g.,  Hibberd  and  Peregrine  (1979).  The  beach  slopes  are  so  gradual 
that  the  waves  break  many  wavelengths  off  shore.  Thus,  in  the  context  of  the 
present  theory,  run-up  occurs  in  a  second  inner  or  near-shore  region. 

In  the  next  section,  we  introduce  the  Korteweg-deVries  equation  and  discuss 
some  of  the  limitations  on  its  use  in  water  wave  problems.  In  Sections  4  and 
5  we  derive  the  main  results  of  the  report  and  compare  these  to  other  studies 
of  this  equation.  In  Section  6,  we  compare  the  results  of  this  asymptotic 
theory  to  numerical  solutions  of  the  Korteweg-deVries  equation. 
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2.  The  Variable  Coefficient  Korteweg-deVries  Equation  and  its  Limitations. 

We  will  consider  one-dimensional  surface  waves  of  length  A  and  amplitude  a 
propagating  into  water  of  variable  depth.  The  undisturbed  depth  of  the  water  is 

H  =  hQh  (5), 

where  x  is  the  distance  measured  parallel  to  the  free  surface  in  the  direction 
of  propagation,  hQ  is  the  depth  at  x  =  0,  i.e.,  h(0)  =  1,  and  L  is  a  scale 
length  measuring  the  rate  of  change  of  the  bottom  topography,  see  Figure  1. 

When  the  parameters 

h 

-  _  a  _  o  ,  .  _  A 

€  =  iT  »  ^  =  T"  and  5  =  l 

o 

are  all  small,  Johnson  (1973a),  Svendsen  and  Hansen  (1978)  and  Cramer  (1980) 
have  all  shown  that  the  wave  motion  is  governed  by  the  following  equation 


1 

2 


,3,2  2  h 

n„  +  t  €  h  nn  +^71) 


,  6  h’  „ 
6  "TXT  4  h  n  °» 


(1) 


where  n  =  hQ  €  q  and  n  is  the  displacement  of  the  fluid  surface  from  its 
undisturbed  position.  Furthermore, 

x  i  -7  ,  h  =  h(x6)  and  h'(xo')  5  7^7 

A  0X0 

This  equation  has  been  written  in  a  coordinate  system  which  moves  at  the  speed 
of  a  linearized  shallow  water  wave;  thus, 

rx& 

ap 

-  t. 


.  1 

T  =  <$  o 


dp 
/h(pT 


where 


t  =  >8ho  7 


and  t  is  the  time  and  g  is  the  acceleration  of  gravity.  In 


the  remainder  of  this  section  we  will  discuss  possible  limitations  on  the  use 
of  this  equation. 

In  using  the  above  equation  it  is  crucial  that  the  ratios  represented  by  €, 
y  and  5  remain  small;  that  is,  the  amplitude  must  remain  small  compared  to  the 
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depth  and  the  length  must  be  large  compared  to  the  depth  and  small  compared  to 
the  scale  length  of  the  bottom  variation.  These  conditions  can  fail  as  the  shore 


is  approached.  However,  we  only  make  use  of  this  equation  up  to  the  point  at 
which  the  wave  breaks.  We  therefore  only  need  to  show  that  these  ratios  are  small 
at  the  point  of  minimum  depth,  i.e.,  the  depth  at  which  breaking  occurs.  To 
provide  a  simple  yet  convincing  argument  that  such  waves  exist,  we  shall  discuss 
the  particular  depth  variation  given  by 


h  =  (1-xS) 

4 

where  p  >  y.  We  first  consider  the  small  amplitude  assumption.  The  ratio  of  the 
amplitude  of  the  wave  to  the  local  depth  is  given  by 

n 


h  h(y) 
o  L 


6n(x,t) 

h(xo) 


In  the  following  sections  we  will  see  that  n  *  0(h  )  at  breaking.  In  order 

for  the  small  amplitude  assumption  to  be  valid  at  breaking  we  must  therefore 


require 


€<<  h. 


5 

4 

b  ’ 


where  h,  is  the  depth  at  which  the  wave  breaks.  It  is  shown  later  that,  for 
b 


the  above  depth  variation, 


hb  =  0(f)7P"4  =  °(1) 


Thus,  at  breaking,  the  wave  amplitude  will  be  small  compared  to  depth  provided 

1  2^£ 

€<<  6  <<  C  ^  ^  ; 

the  first  condition  will  be  seen  to  be  necessary  in  order  to  induce  breaking 
rather  than  fissioning.  The  above  order  relation  is  self-consistant  provided 
2  2-p  4 

—  — ^  <  1,  i.e.,  p  >  ■=■  ;  this  is  seen  to  come  naturallv  out  of  the  analvsis 
op/’ 

later.  Note  that  the  amplitude  to  depth  ratio  is  of  order 


2  £^2  5p 

(€5  P  6)7p'4 


at  breaking.  Thus,  when  p  >  2,  this  ratio  is  always  small  and  we  only  need 
to  require 

€«  6  «  1  . 

We  now  consider  our  long  wave  assumption.  As  the  wave  approaches  shore, 
the  linearized  shallow  water  theory  predicts  that  the  length  of  the  wave 
decreases.  However,  the  rate  at  which  the  depth  decreases  is  much  faster  than 
this  shortening  of  the  wave.  Until  the  wave  comes  very  close  to  shore,  the 
long  wave  approximation  is  actually  getting  better  rather  than  worse.  The 
main  effect  may  be  illustrated  using  the  linearized  shallow  water  theory.  If 
we  ignore  reflected  waves  (these  will  only  increase  the  length  of  the  wave, 
thereby  making  our  approximation  better),  this  theory  may  be  used  to  show  that 
the  ratio  of  the  depth  to  the  length  of  the  wave  is  given  by 


izR 

2p 


provided 

5  <<  h 

In  order  to  ensure  that  this  long  wave  assumption  is  not  violated,  we  therefore 
require 


6  << 


"b 


IzR 

2p 


Because 


bb 


4P  .. 

ocf) 7p-4 


we  have 


2  2^2 


€«  5  <<  €  5  P  , 

which  is  the  same  as  that  obtained  for  the  small  amplitude  assumption. 

We  now  consider  the  assumption  of  a  gradually  varying  bathymetry.  If  we 
again  use  the  linearized  shallow  water  theory  and  neglect  reflections,  we  find 
that  the  ratio  of  the  length  of  the  wave  to  the  bathymetry  length  scale, 


5 


p 


-1 

1 1  dhi 

I'jJ'  2=1  ,  is  given  by 


provided  <5  <<  h 


ir£ 

2p 


Thus , 


6 


h  2p 


the  conditions  imposed  earlier,  viz. 


_2=P 

s  «  S  2p 


I 

i.e.,  €«  5  «€5  p  , 

are  sufficient  to  ensure  that  this  ratio  remains  small  even  at  the  breaking  point. 

Thus,  we  see  that  the  wave  will  always  break  before  the  small  amplitude, 

long  wave  and  gradual  depth  variation  assumptions  fail  provided  we  require 

2  2-£ 

€  «  6  «  €5  p  «  1. 

If  p  >  2  these  assumptions  will  always  be  valid  at  breaking  and  we  do  not  need 
to  impose  any  conditions  other  than  those  of  Sections  3-5,  i.e., 

€<<<$<<  1. 

It  should  be  noted  that  these  results  were  obtained  for  a  special  form  of  the 

depth  variation;  it  is  expected  that  this  discussion  is  easily  extended  to 

more  general  bathyraetries.  The  conditions  given  here  are  related  to  the 

local  conditions  given  by  Varley,  Venkataraman  and  Cumberbatch  (1971) . 

Equation  (1)  also  assumes  that  we  may  neglect  the  effect  of  reflected 

waves.  The  above  discussion  shows  that  the  length  ->f  the  wave  remains  small 

compared  with  the  length  scale  associated  with  the  depth  variation.  It  is 

therefore  clear  that  the  reflected  waves  generated  by  the  depth  variation  will 

be  small,  at  least  for  the  time  period  considered  here.  For  wavetrains,  we 

expect  the  effect  of  reflections  to  be  diminished  even  further  by  destructive 

interference,  see,  e.g.,  Kajiura  (1961).  Previous  numerical  and  experimental 

studies  also  provide  support  for  this  conclusion.  Peregrine  (1967)  has  used 
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the  Boussinesq  equations  to  describe  a  solitary  wave  as  it  propagates  over  a 
gently  sloping  beach;  the  resultant  reflected  wave  was  observed  to  be  at  least 
an  order  of  magnitude  smaller  than  the  incident  wave.  Waves  which  break  by 
plunging  without  significant  reflection  have  also  been  observed  in  the  laboratory 
studies  of  Street  and  Camfield  (1966)  and  Galvin  (1968). 

We  conclude  that  the  Korteweg-deVries  equation  (1)  will  provide  a  valid 
model  for  shoaling  water  waves  under  relatively  mild  restrictions  on  the  initial 
wave  shape  and  bottom  topography.  It  should  be  noted  that  the  theory  here  only 
attempts  to  calculate  the  conditions  associated  with  the  breaking  of  the  wave. 

In  the  final  stages  of  the  approach  to  shore,  most  of  the  assumptions  made  here 
fail  and  one  must  resort  to  more  complicated  theories  capable  of  treating  the 
problems  of  run-up  and  backwash. 
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3.  Problem  Formulation 


In  this  report  we  will  obtain  approximate  solutions  to  equation  (1)  subjec 
to  the  initial  condition 


n(0,x)  =  F(x)  . 


If  we  define 


x  =  x<$  and  u(x,x) 
the  initial  value  problem  reads 


(x)rj(x. 


t) 


u_ 

X 


+  £  h2U 
6  TXT 


+  ^sh 


l 

4 

uu 

x 


0 


(2) 


u(0,x)  =  F(x) , 

2 

where  h  =  h(xo)  =  h(x) ,  v  5  7—  and  e  i  t. 

o  o 

On  Figure  1  we  have  sketched  typical  bottom  topographies.  Figure  la 
depicts  a  closed  beach;  here  h(l)  =  0.  In  Figure  lb  a  shelf  is  sketched; 
here  h(x)  =  d  for  x  _>  1.  In  the  numerical  studies  described  below  we  will 
always  take  h  =  1  for  x  <.  0.  It  should  be  noted  that  the  length  scales  have 
been  distorted  for  the  purposes  of  illustration. 

In  this  report  we  will  consider  the  case  where  e  is  relatively  small, 


x.e. , 


€<<  5  <<  1. 


The  parameter  v  is  taken  to  be  of  order  one,  i.e.,  at  most  bounded.  Inspec¬ 
tion  of  equations  (2)  suggests  that  the  wave  motion  may  be  described  by 
the  linearized  version  of  this  equation 


1 

L  v  ,2 
u_  +  t*  h  u 
x  6 

u(0,x)  = 


XXX 

F(x) 


=  0 


at  least  for  X  i  1,  i.e.,  h  i  0.  Of  course,  if  v  is  also  small,  e.g., 
2 

u  =  0(s),  then  the  wave  motion  is  described  by  the  classical  Green's 
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law  formula  u(x,t)  :  F(t),  see,  e.g..  Lamb  (1932),  However,  as  the 
shore  is  approached,  that  is,  as  h  +  0,  the  nonlinear  term  is  expected 
to  become  larger  and  larger  and  the  linear  solution  is  therefore  expected 
to  break  down  near  the  shore.  An  inspection  of  the  dispersive  term  in 
(2)  suggests  that  it  will  become  smaller  and  smaller  as  h  -*■  0.  We  would 
therefore  expect  that  the  approximation  to  (2)  which  is  valid  near  the 
shore  is 

l 

3  "  4  _  A 

u^  +  meh  uu  ,  0. 
x  2  x 

Of  course,  the  initial  condition  for  this  equation  is  simply  the  incoming 
linear  solution.  This  equation  is  immediately  recognized  as  a  version  of 
the  inviscid  Burger  equation  or  Airy's  equation.  Solutions  to  this  equa¬ 
tion  steepen  at  the  front  face  of  the  wave  and  eventually  become  triple¬ 
valued.  Although  the  theory  must  fail  where  the  surface  slope  becomes 
infinite,  the  usual  interpretation  of  this  phenomena  is  that  the  actual  wave 
will  be  observed  to  topple  over  and  break  in  a  similar  fashion,  see,  e.g., 
Stoker  (1957)  .  These  remarks  are  expected  to  hold  not  only  for  a  closed 
beach  but  for  a  shelf,  see  Figure  lb,  provided  d  is  sufficiently  small. 

This  heuristic  description  of  the  wave  evolution  is  shown  to  be  an 

accurate  one  in  the  following  sections.  In  order  to  illustrate  the  basic 

theory,  the  case  of  v  *  0(c)  =  o(l)  is  described  in  detail  in  Section  4;  this 

will  also  enable  us  to  quickly  extract  results  appropriate  for  solitary 

2 

waves.  Because  €  and  y  are  small  compared  to  5  it  will  be  referred  to 
as  weakly  dispersive  and  weakly  nonlinear.  In  Section  5,  shorter  waves 
having  v  *  0(1)  are  discussed.  Because  these  exhibit  some  dispersion 
before  entering  the  near-shore  region  and  breaking,  this  will  be  referred 
to  as  the  dispersive  case. 


9 


t 


4.  Weakly  Dispersive  and  Weakly  Nonlinear  Waves 

In  this  section  it  will  be  assumed  that  both  e  and  v  are  small.  We 
assume  an  outer  or  off-shore  expansion  of  the  form 

u  *  UQ  +  0(e,v) , 

which  is  assumed  to  be  valid  for  all  x  =  0(1),  x  i  1.  If  we  substitute  this 
in  equation  (2)  and  retain  only  lowest  order  terms,  we  find  that 

Uq(x,t)  =  F(t) , 

which  implies, 

u  =  F(t)  +  0(s, v) , 

for  is  )i  1.  This  can  be  recognized  as  Green's  law. 

An  examination  of  higher  order  terms  shows  that  this  straightforward 
expansion  scheme  breaks  down  in  the  vicinity  of  the  shore,  i.e.,  x  =  1.  We 
therefore  need  to  find  a  second  approximation  to  the  Korteweg-deVries  equation 
in  (2)  which  is  valid  near  the  shore.  To  do  this  we  define  a  new  independent 
variable 

* 5  ¥  >  (4) 
where  A  =  o(l)  is  a  scale  factor  to  be  determined  later.  Note  that  &  <  0  when 

x  <  1.  In  the  usual  way,  we  will  assume  that  &  =  0(1)  in  the  inner  or  near-shore 

region;  the  scaling  1  essentially  gives  the  nearness  of  x  to  1.  We  will  also 

need  to  write  out  explicitly  the  expansion  of  h  in  the  vicinity  of  x  *  1.  This 

will  be  written 

h(x)  =  h(l+Aft)  ~  AahU)  +  o(Aa)  (5) 

for  small  A  and  $  =  0(1) .  The  exponent  a  >  0  and  the  form  of  the  function 
h(k)  will  depend  on  the  particular  behavior  of  h  as  x-*-l.  For  example,  for  a 
closed  beach  and 

h  -  (l-x)n  as  k-^-l, 
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we  have  a  =  n  and  h(x)  =  (-x)  .  For  che  shelf  given  by 

1  for  x  <  0 

1+d  ,  1-d 


h(x)  = 


2  +  — cosllx 


where  d  is  assumed  small,  we  have  a  =  2  and 


for  0  <  x  <  1 
for  x  >  1 


h(x)  = 


*  j[2it2 

d  +  ■—  for  St  ^  0 
4 


'  d  for  &  >  0 

where  d  i  A  2d.  The  most  interesting  case  is  when  d  =  0(1) . 

If  we  now  use  (4)  to  change  variables  in  (2),  substitute  the  expansion 
(5)  and  the  following  inner  or  near-shore  expansion 


u  =  u  +  o(l) 
o 


in  (2),  we  then  have 


i ,  a  1  7-7 

1+  •?  1-  t  cu  ~r 

u+vA  hu  +  -  ei  huu=  o(l). 

o.  ottt  2  o  0T 

X 

For  v  =  o(l) ,  e  =  o(l)  and  a  >  0,  the  only  reasonable  balance  of  terms  is 
obtained  when  , 


7a-4  4 

A  =  e  and  a  >  y. 

In  the  former,  the  order  symbols  have  been  dropped  for  convenience.  The 

Korteweg-deVries  equation  therefore  becomes 

-7 
3  -  4 

u  ,  +  y  h  uu  =0, 
ox  2  o  ox 

to  lowest  order.  The  general  solution  to  this  equation  is  easily  obtained 
by  the  method  of  characteristics;  this  may  be  written  in  the  following 
parametric  form 

u  U,t)  =  G (x  ) 

O  03 
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and 


T  =  T»  +  1  G(tJI(x), 

where  I(x)  =  Hi  (p)dp  .  (6) 

—  00 

Here  G  is  an  unknown  function  to  be  determined  by  matching  to  the  linear 

solution.  Alternatively,  one  could  recast  this  as 

u  (x,  t)  =  G(x-  4  GI (ft) )  , 
o  Z 

which  just  follows  from  substituting  from  the  second  parametric  equation 
into  the  first. 

Our  description  of  the  near-shore  solution  will  be  complete  once  G  is 
specified.  This  will  be  obtained  by  matching  the  near-shore  solution  to 
the  off-shore  solution  through  the  formalism  of  the  method  of  matched 
asymptotic  expansions,  see,  e.g.,  Van  Dyke  (1975).  We  first  consider  the 
near-shore  solution 

u  =  G(r  —  |  GI (ft) )  +  o ( 1 )  . 

When  this  is  recast  in  the  off-shore  variables  and  expanded  for  e-»0  and 
x  fixed,  we  have 

u  '  G(t)  +  o(l)  , 

_4 _ 

X- 1  7  qi  —  A 

because  I  (ft)  =  I(— ^)  -*  0  as  A  =  e  -*■  0.  The  method  of  matched  asymptotic 
expansions  then  requires  that  this  match  the  off-shore  solution 

u  =  F(t)  +  o(l)  . 

As  a  result  G(p)iF(p),  for  all  p,  and  the  near-shore  solution  therefore  becomes 

u  -  F(r  -  |  FI (ft) )  +  o(l)  , 
or,  in  terms  of  the  characteristic  lines, 

U (ft , T )  *  F(I J  +  0(1) 
on  r  -  t  +  4  F(t  )I(ft)  . 

CO  X 
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It  is  well-known  that  some  of  the  characteristics  will  converge  on  themselves 
leading  to  a  steepening  of  the  wave  at  the  front  face.  The  break  point 
will  be  defined  as  the  point  at  which  the  slope  of  the  wave  first  becomes 
infinite.  The  equation  of  the  characteristic  lines  (7)  may  be  used  to  show 
that  this  is  given  by  the  following  implicit  formula: 

=  3  IF'  j  ’  ( 

1  I  max 

dF 

where  F'(^)  =  ~  .  Here  we  note  that  higher  order  matchings  have  been  carried 
out;  these  produce  no  changes  in  the  basic  lowest  order  results. 

In  summary,  we  have  derived  approximate  solutions  to  the  initial  value 
problem  (2)  for  small  v  and  small  e.  We  have  found  that  for  all  x  i  1  the 

behavior  of  the  wave  is  given  by  Green's  law,  i.e., 

u  =  F(t)  . 

Close  to  the  shore,  this  description  is  no  longer  valid  and  it  is  seen  that 
nonlinear  effects  play  a  dominate  role.  The  resultant  solution  and  formula 
for  the  break  point  are  given  by  equations  (7)  and  (8).  The  relation  of 

the  near-shore  variable  k  to  the  variables  x  and  x  are 


7a-4  „ 

x=xo=l+Ax=l+c  x. 


where  a  >  y  .  Thus,  breaking  occurs  when  x-1  =  0(e  )  and  the  depth  hj. 


at  which  breaking  occurs  is 


hb  -  0(e?“-4), 


where  (5)  has  been  used.  At  breaking,  the  actual  wave  height. 


n  =  €hQh  u, 


.  -  ,  ,  -7a-4  ,7a-4 

is  theretore  ot  order  he  o 

o 

It  is  interesting  to  note  that  dispersive  effects  play  no  role  in  the 


wave  evolution,  at  least  in  the  case  considered  here.  In  the  off-shore  or 
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outer  region  they  are  negligible  for  the  same  reason  nonlinear  effects  are; 
that  is,  they  are  only  important  over  distances  which  are  large  compared  to 
the  distance  to  shore.  Near  the  shore,  h  is  small  and  the  dispersive  term  in 
(2)  is  negligible  compared  to  the  nonlinear  term.  In  fact,  we  see  that  this 
case  may  be  modeled  by  the  nondispersive  version  of  (2),  viz., 

I 

,  3  . “  4 

u_  +  —  eh  uu  =0 
x  2  t 

u(0,x)  =  F(t). 

To  lowest  order,  the  solutions  to  this  problem  are  identical  to  those 
obtained  here. 

It  should  also  be  noted  that  the  results  obtained  here  are  valid  for 
solitary  waves,  i.e.,  €  =  0(u2) .  Because  €<<&<<  1  we  might  refer  to  these 
as  strongly  shoaling  solitary  waves.  In  fact,  in  terms  of  the  ordering  of 
the  parameters  €,  u2  and  5,  this  is  exactly  the  case  studied  by  Tappert 
and  Zabusky  (1971),  These  authors  have  considered  the  case  of  a  shelving 
solitary  wave  and  related  the  number  of  solitons  generated  to  the  depth 
ratio  of  the  shelf.  The  problem  discussed  in  this  section  may  be  applied 
to  shelving  waves  and  it  is  clear  that  for  sufficiently  small  shelf  depths 

the  wave  will  break  without  fissioning.  In  the  example  discussed  earlier 

4_ 

€  5 

in  this  section  the  shelf  depth  must  be  of  order  (-j)  in  order  for  breaking  to 
occur.  Thus,  this  theory  provides  limits  on  that  of  Tappert  and  Zabusky  and,  in 
particular,  limits  the  number  of  solitons  which  can  be  generated  by  the  shelving 
process . 

In  Section  6,  the  results  obtained  here  are  compared  to  numerical  solutions 
of  (2) .  It  is  also  desirable  to  compare  these  results  with  calculations  which 
include  reflected  waves.  One  such  study  is  due  to  Peregrine  (1967).  The 
appropriate  Boussinesq  equations  were  derived  and  solved  numerically.  The 
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resultant  evolution  was  interpreted  by  Peregrine  to  include  breaking  near 
the  shore.  The  present  study  neglects  reflected  waves;  the  fact  that  these 
remain  small  up  to  breaking  is  certainly  suggested  by  Peregrine's  calculations. 
Furthermore,  the  wave  growth  appears  to  be  in  reasonable  agreement  as  well. 

To  see  how  this  fits  in  with  previous  theories,  consider  a  solitary  wave 
of  fixed  amplitude  and  length,  i.e.,  fixed  €  and  u,  €=  0(u2).  If  we  then 
let  this  approach  a  very  gradual  beach,  e.g.,  €  »  5  or  e  >>  1,  we  obtain 
the  gradually  varying  solitary  wave  solution  of  Johnson  (1973b).  If  we  let  the 
same  wave  approach  a  second  beach  which  is  steeper,  say  &  =  0(()  or  i  =  0(1), 
we  find  that  the  dispersion,  modulation  and  nonlinearity  all  act  with  equal 
strength  and  the  wave  fissions  as  described  by  Madsen  and  Mei  (1969).  For  an 
even  steeper  beach,  i.e.,  6  >>  €  or  e  =  o(l),  the  wave  will  come  relatively 
close  to  shore  without  significant  distortion.  Once  the  wave  is  very  near 
shore  the  dispersion  is  very  weak  and  the  wave  must  break  as  described  by  the 
present  theory. 
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We  now  consider  shorter  waves  which  exhibit  a  dispersive  behavior  before 


breaking  occurs.  The  solution  to  the  initial  value  problem  (2)  will  now  be 

outlined  for  the  case  v  =  0(1)  and  e  *  o(l) ,  that  is,  u2  *  0(6)  and  €=  o(5) . 

We  now  assume  an  off-shore  expansion 

u  =  U  +  o  ( 1) , 
o 

where  there  should  be  no  confusion  between  the  notation  in  this  section  and 
the  previous.  The  initial  value  problem  (2)  therefore  becomes 


U  ,  +  r  h *(x)U  =  0 

OX  6  OTTT 


Uo(0,t)  =  F(t), 

to  lowest  order.  We  can  transform  this  into  a  problem  in  wate*-  of  constant 
depth  by  the  following  transformation 


5 


Z  = 


h^(p)dp. 


which  yields 


0,+j0 

OZ  6  OTTT 


Uq(0,t)  »  F(t). 


(9) 


General  solutions  to  (9)  may  be  obtained  by  transform  techniques.  Because 
v  =  0(1)  these  are  strongly  dispersive  solutions.  If  we  assume  that 

F(t)  =  cos  x, 

an  exact  solution  to  (9)  is 

Uq(x,t)  =  cos  T, 


where  +  J  hH(p)dp. 


■  •*«*<*■ 
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As  in  Section  4,  we  will  need  a  second  approximation  to  (2)  near  the 
shore.  By  the  same  procedure  as  used  in  Section  4  we  find  that 

uo(*.t)  =  G(tJ 

on  t  =  rx  +  |  G(TjI(ft), 

x  _  7 

where  I(&)  =  h  ^(p)dp 

—00 

and  G  is  an  unknown  function  to  be  determined  by  the  matching.  Here  the 
near-shore  expansions 

u  =  u  +  o(l) 
o 

4 

4-7a 

*  =  (x-1)  e  *  0(1) 

have  been  used.  Again,  one  should  not  confuse  the  notation  of  this  section 
and  the  last. 

Because  of  the  dispersion  present  in  the  off-shore  solution,  see,  e.g., 
(9),  it  is  clear  that  the  result  of  the  matching  will  not  yield  G  =  F. 
Essentially  G  will  be  the  solution  to  (9)  evaluated  at  x=l.  For  example, 
when  F(t)  =  cos  t  the  near-shore  solution  is  found  to  be 

u  (x,t)  =  cos  T 

O  00 

on  3 

T  =  T»  +  2  I(8)  C0S  T»’ 

v  f  u 

where  +  -g  h^(p)dp,  or,  equivalently, 

o 

u  (11, t)  *  cos(x  +  h^(p)dp  -  ^  1(A)  cos  T  ) 

O  0  Z  50 

o 

3 

on  t  =  t  +  •=■  I(k)  cos  T  . 

»  2.  OO 
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A  breaking  condition  may  also  be  derived;  in  this  case,  and  will 
depend  on  v  through  T  .  As  in  Section  4,  higher  order  matchings  have  been 
carried  out;  these  show  that  the  solutions  presented  above  are,  in  fact, 
the  lowest  order  approximations  to  (2)  in  the  appropriate  regions. 

Thus,  we  see  that,  although  the  dispersion  plays  a  major  role  in  the 
off-shore,  x  =  0(1),  evolution  of  the  wave,  it  is  negligible  near  the  shore. 
The  breaking  point  depends  on  v  or  u2  only  insofar  as  the  "incoming  wave" 
does.  It  is  clear  that  the  results  sketched  here  will  contain  those  of 
Section  4  when  v  is  small. 


t 


6.  Numerical  Results 


In  this  section  we  compare  numerical  solutions  of  (2)  to  the  asymptotic 
solutions  derived  in  this  report.  We  shall  concentrate  on  the  behavior  of  solitary 
waves  and  therefore  the  results  of  Section  4  will  be  of  interest.  The  finite 
difference  scheme  developed  by  Zabusky  and  Kruskal  (1965),  Vliegenthart  (1971) 
and  Johnson  (1972)  was  applied  to  (2)  for  various  values  of  v  and  z.  As  in 
the  previous  investigations,  the  differential  equation  was  replaced  by  the 
difference  formula 

v  ,  %  Ax 


ui,j+l  =  Ui,j-1  6  h  (At!3  (ui+2,j  "  2ui+l,j  +  2ui-l,j  "  Ui-2,j) 
7 


■f  h.  ^  ~  (u.j,  .  +  u.  .  +  u  .)(u  .  .  -  u  .) 

2  2  At  i+l, j  1,3  l-l, j  i+l,j  1-1,3 


where  i  ai.d  j  correspond  to  the  t  and  x  directions,  respectively,  At  and  Ax 


are  the  step  sizes  in  the  t  and  x  directions  and  tu  =  h(jAx).  The  first  step 


was  calculated  by  a  simple  uncentered  scheme.  For  this  problem,  the  stability 
criterion  reads 


Ax 

At 


l 

4 


% 


2  e 


,  h. 

+  1  J  , 
3  TatT7 


v)  <  1 


It  was  found  that  step  sizes  of  At  =  0.1  and  Ax  =  0.0001  satisfy  the  above 
criterion  and  provided  satisfactory  results  in  all  the  calculations  described 
below. 

We  can  expect  differences  between  the  asymptotic  results  of  Section  4  and 
those  obtained  by  the  finite  difference  scheme  for  two  reasons.  The  first 
is  that  the  results  of  Section  4  are  only  valid  in  the  limit  of  small  e.  As 
we  would  expect,  the  agreement  between  the  two  solutions  becomes  better  and 
better  as  e  approaches  zero.  The  second  source  of  discrepancy  is  due  to  the 
numerical  dispersion  inherent  in  the  finite  difference  scheme.  For  the  sake 
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Em 


of  illuscration,  consider  Che  following  initial  value  problem  involving 

che  constant  coefficient  Burger's  equation. 

v  +  w  =0 
x  T 

v(0,t)  =  F(x) 

Exact  solutions  to  this  problem  are  easily  obtained  by  the  method  of  characteristics. 
It  is  well-known  that  the  waveform  steepens  until  the  slope  at  some  point  on  the 
front  face  of  the  wave  becomes  infinite.  After  this  time,  the  solution  is  seen 
to  become  triple-valued.  However,  when  this  problem  is  solved  using  the  finite 
difference  scheme,  the  predicted  behavior  is  quite  different.  Although 
during  the  initial  stages  the  behavior  is  similar  to  the  exact  solution  and 
although  the  wave  steepens  noticeably,  the  wave  slope  never  becomes  infinite. 

In  fact,  the  wave  is  observed  to  fission  in  a  manner  similar  to  that  of  the 
constant  coefficient  Korteweg-deVries  equation  and  is  seen  to  emit  a  large 
number  of  very  narrow  spike-like  solitons.  If,  in  a  Korteweg-deVries  equation, 
the  physical  dispersion  is  weak  enough,  e.g.,  of  the  order  of  the  numerical 
dispersion,  erroneous  results  are  also  expected.  Thus,  near  breaking,  some 
care  should  be  taken  when  comparing  the  results  of  our  asymptotic  theory  to 
the  numerical  results.  A  crude  measure  of  this  error  was  obtained  by  comparing 
exact  and  finite  difference  solutions  to  the  appropriate  Burger's  equation;  this 
is  discussed  later. 

The  waves  considered  were  of  the  solitary  wave  type.  The  initial  condition 
used  was 

F(t)  h  sech2  t  . 

In  order  to  iLlustrate  the  strength  of  the  shoaling  we  have  chosen 

4 

v  -  3  E  . 

As  a  result,  this  wave  will  fission  in  water  of  uniform  depth.  However,  we 
will  see  that  the  strength  of  the  shoaling  is  such  that,  even  in  this  case. 
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the  fissioning  cannot  occur  and  the  wave  breaks  as  described  in  Section  4.  For 


the  sake  of  simplicity,  the  beach  slope  was  taken  to  be  constant.  That  is, 

!1  x  <  0 

1-x  0  ^  x  <  1; 

thus,  the  shore  is  located  at  x  =  1  or  x  =  6  \ 

In  each  case,  the  asymptotic  solution  has  been  calculated  from  equations 
(6)  and  (7) .  For  the  beach  and  wave  profile  described  above,  the  solution 
may  be  calculated  from 

u(k,x)  =  sech2  t 

00 

_  2  do) 

x  =  t  +  2(-£)  ^  sech2  x  , 

00  00 

4 

3 

where  x^  may  be  regarded  as  a  parameter  and  £  =  (x-l)e  .  The  breaking  point 
is  calculated  from  (8)  to  be 

4 

\  i  1-x,  =  (2e|F*  |  )3 

b  b  11  max 


or 


hK  =  (- 


3/3~ 


t)3. 


since  |Ff 


‘max 


3  JT 


for  this  case. 


In  Figure  2,  the  evolution  of  a  wave  having  e  *  0.1  has  been  depicted.  The 
solid  lines  are  the  finite-difference  solutions  and  the  dotted  lines  the  asymptotic 
theory.  Figure  2(a)  is  the  waveform  at  S  =  0,  Figure  2(b)  is  the  waveform  at 
k  =  0.9  and  Figure  2(c)  is  the  waveform  at  breaking  x^  *  0.918.  The  discrepancy 
is  expected  to  be  due  to  the  small  e  approximation.  Figure  3  depicts  the 
evolution  for  a  smaller  value  of  z;  t  *  0.05.  As  we  would  expect,  the  agreement 
is  much  better;  again,  the  error  is  within  the  bounds  of  the  small  e  approximation. 
As  predicted  by  the  theory  of  Section  4,  there  is  very  little  steepening  of  the 


I 


wave  until  very  near  the  shore,  see,  e.g.  Figure  3(b).  In  Figure  4,  the  waveform 
corresponding  to  e  =  0.025  is  shown.  Here  the  agreement  is  quite  good,  even  at 
the  breaking  point;  see  Figure  4(c).  The  main  source  of  discrepancy  is  expected 
to  be  due  to  numerical  dispersion.  As  a  check  on  some  of  the  assertions  made  in 
this  report,  we  have  recalculated  the  waveform  using  the  Burger's  equation: 

1_ 

.  3  '  4  . 

u_  +  -r  eh  uu  =0 

X  2  T  (; 

u(0,t)  =  sech2T, 

where  e  =  0.025.  We  first  obtained  the  exact  solution  to  this  problem  by  the 
method  of  characteristics  and  then  an  approximate  solution  was  calculated  using 
the  above  finite  difference  scheme.  In  Figure  5(a) ,  we  have  compared  the 
asymptotic  solution  (10)  to  the  exact  solution  of  (11);  as  asserted  at  the  end 
of  Section  4,  we  see  that  these  agree  very  well.  In  Figure  5(b),  we  have 
compared  the  finite  difference  solutions  of  (2)  and  (11).  These  are  also  in 
agreement  suggesting  that  the  physical  dispersion  in  (2)  is  very  small;  at 
least  in  the  problem  considered  here.  Because  the  only  difference  between  the 
finite  difference  and  exact  solutions  to  (11)  is  due  to  numerical  dispersion, 
a  comparison  of  all  four  plots  in  Figure  5  strongly  suggests  that  the  discrepancy 
in  Figure  4(c)  is  primarily  due  to  the  numerical  dispersion  inherent  in  the 
finite  difference  scheme. 


7.  Conclusion 


This  report  describes  a  particular  class  of  solutions  to  the  variable 
coefficient  Korteweg-deVries  equation.  For  the  most  part,  nonlinear  effects 
were  negligible  and  the  wave  evolution  was  described  by  a  linear  and,  in  some 
cases,  dispersive,  equation.  Because  of  the  amplification  and  effective 
lengthening  of  the  wave,  we  find  that  near  the  shore  dispersive  effects  are 
negligible  and  the  nonlinear  steepening  effects  predominate;  this  leads  to  a 
breaking  of  the  wave.  A  careful  consideration  of  higher  order  terms  has  shown 
that  the  perturbation  and  matching  procedure  is,  in  fact,  correct;  furthermore, 
the  numerical  solutions  presented  here  give  further  support  for  this.  Thus, 
we  expect  that  when  exact  solutions  to  (2)  are  obtained,  these  will  approach 
ours  in  the  limit  of  vanishing  e  H  ~  . 

O 

We  have  also  provided  estimates  for  the  break  point  and  the  wave  height 
at  breaking.  In  the  case  of  solitary  waves,  or,  more  generally,  £=  0(y2),  we 
find  that  these  depend  on  the  parameters  €,  u  and  5  only  through  e  =  j  .  For  the 
general  dispersive  case,  i.e.,  v  =  0(1),  there  is  an  additional  dependence  on 

u2 

v  =  —  through  the  initial  condition  for  the  near-shore  problem. 

It  is  expected  that  the  results  obtained  here  will  also  provide  a  good 
approximation  to  actual  waves  on  a  beach.  While  we  can  find  no  numerical  or 
experimental  studies  which  permit  direct  comparison,  we  find  that  our  results 
are  in  qualitative  agreement  with  the  numerical  study  of  Peregrine  (1967)  and 
the  laboratory  studies  of  Street  and  Camfield  (1966)  and  Galvin  (1968);  in 
particular,  the  results  of  Section  4  correspond  to  the  plunging  breaker  regime 
of  Street  and  Camfield. 

This  work  was  initiated  and  carried  out  while  the  first  author,  M.S.C., 
participated  in  the  ASEE-NAVY  Summer  Faculty  Research  Program  at  the  NRL.  The 
second  author  B.M.  was  supported  by  NRL  under  contract  number  N00173-80-M-7172. 
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I 


Figure  1 


1 


(a) 


(b) 

Typical  Sathymetrles .  The  undisturbed  water  level  is  taken  to  be 
y  ■  0  and  L  may  be  considered  to  be  the  length  scale  characteristic  of 
the  depth  variation,  (a)  Closed  beach,  (b)  Shelf;  d  gives  the  ratio 
of  the  shelf  depth  to  the  deep  water  depth. 
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Wave  Evolution  for  Solitary  Wave  Protile  u(0,r)  ■  sechz  r  and  ■ 
v  *  0.133.  Asymptotic  theory  -  -  -  -  ,  Finite  Difference  Solution 
(a)  *  -  0.  (b)  X  -  0.9,  (c)  S  -  0.918. 


(a) 


T 


r 

(c) 


Figure  3:  Wave  Evolution  for  e  «  0.05,  v  ■  .0667.  Asymptotic  Theory  -  -  -  -  , 

Finite  Difference  Solution  -  .  (a)  x  »  0,  (b)  *  *  0.9,  (c)  X  «  .967. 
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(b) 

Figure  5:  Comparison  of  Various  Theories  for  u(C,t)  «  sech2  t,  e  »  0.025, 

v  ■  .0333,  *  ■  .987.  (a)  Asymptotic  Solution  -  ,  Exact  Solution 

of  Burger's  Equation  -  -  -  -  ,  (b)  Finite  Difference  Solutions: 

Full  Korteweg-deVries  Equation  - -  ,  Burger's  Equation  -  -  -  . 
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